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This  paper  is  concerned  with  the  modelling  of  the  return  to  isotropy  part  of 
the  pressure  strain  term  in  homogeneous  anisotropic  turbulent  flows. 
Analytical  solutions  of  the  transport  equations  of  the  invariants  of  the 
anisotropy  tensor  as  well  as  that  of  turbulent  kinetic  energy  as  a  function  of 
the  natural  time  of  decay  are  provided  and  discussed.  Principal  components  of 
the  Reynolds  stresses  are  obtained  from  the  solution  of  a  cubic  equation  which 
involves  the  invariants.  It  is  shown  that  current  models  based  on  Rotta's 
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b.  STATEMENT  OF  WORK 

The  objective  of  this  research  is  to  develop  a  more 
general  and  universally  applicable  model  of  the 
pressure-strain  correlation  of  the  Reynolds  stress 
transport  equations  with  special  emphasis  on  turbulent 


swi r 1 i ng  f 1 ows .  A 
turbulent  flows 
establishing  the 
constants  on  the 
anisotropy  tensor, 
applied  to  the  case 


systematic  investigation  of  homogeneous 
will  be  continued  with  the  aim  of 
functional  dependence  of  the  model 
second  and  third  invariants  of  the 
The  model  developed  will  first  be 
of  homogeneous  rotating  turbulence  and 


then  to  more  complicated  compressible  flows  and  swirling 
flows  of  practical  importance.  An  experimental  program 
also  is  proposed  to  generate  detailed  measurements  of  both 
the  mean  and  turbulence  quantities  for  model  validation 
purposes . 


c  . 


STATUS  OF  THE  RESEARCH  EFFORT 


It  is  generally  recognised  that  accurate  modelling  of 
complex  practical  flows  would  require  the  use  of  second 
order  closure,  and  that  the  week  point  of  the  closure  is 
the  lack  of  universality  of  the  model  for  the 
p r e s s u r e - s t r a i n  correlation  +  jj  in  the  Reynolds  stress 
transport  equations.  Chou  (1)  showed  that  with  the  aid  of 
a  Poisson  equation  for  the  fluctuating  pressure,  p,  the 
pressure-strain  terms  can  be  re-expressed  in  terms  of  an 
integral  without  the  explicit  appearance  of  p.  The  form  of 
the  integral  suggests  that  three  distinct  processes 
contribute  to  the  pressure-strain  terms  (1,2).  The  first, 
the  'return  to  isotropy'  terra,  +  j  j  i  ,  involves  only  the 
fluctuating  velocities.  The  second,  the  'rapid  term’  or 
the  '  isotropizat  ion  of  the  production  terra'  t i  j  ,  2  ■ 
involves  the  interaction  between  the  mean  strain  and 
fluctuating  velocities.  The  third,  the  wall  proximity 

term.  +ij,w-  expresses  the  effect  of  walls  on  the 
pressure-strain  terra.  The  net  effect  of  these  processes  is 
to  redistribute  the  energy  among  the  Reynolds  stress 
i  components.  but  nothing  can  be  deduced  from  the  Reynolds 

stress  transport  equation  about  the  direction  and  the  speed 
of  the  energy  transfer. 

During  the  course  of  this  work.  substantial 
achievements  have  been  obtained  in  the  modelling  of 
homogeneous  anisotropic  turbulence  without  mean  velocity 
gradient.  Experiments  (3-6)  showed  that  this  kind  of  flow 
embodies  the  general  dynamics  of  the  decay  of  the  energy 
containing  eddies,  a  process  described  by  the  pressure 
strain  tensor.  The  flow  also  enables  tfie  analysis  of 
energy  transfer  and  functional  dependence  of  contants  to  be 
simplified  as  only  the  return  to  isotropy  part  of  the 
pressure-strain  tensor  remains.  Under  these  conditions,  a 
general  form  of  the  current  pressure-strain  model 
incorporating  non-linear  terms  (7-9)  is  as  follows: 


+i  j  = 
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where  bjj  is  the  anisotropy  tensor,  k  the  mean  turbulent 
kinetic  energy,  II  the  second  invariant  of  the  anisotropy 
tensor,  €  the  viscous  dissipation  rate,  djj  the  kronecker 
delta,  ci  and  X  may  be  either  constants  or  functions  of 
some  relevant  parameters,  and  angled  brackets  denote 
averages  over  time.  The  first  term  on  the  right  hand  side 
of  Eq.  (1)  is  the  linear  model  proposed  by  Rotta  (10),  and 
the  second  term  expresses  the  non-linear  effects. 


Closure  of  the  Reynolds  stress  equation  requires 
further  modelling  assumptions  in  its  viscous  dissipation 
term  and  in  constructing  a  length  scale  or  the  dissipation 
transport  equation.  Here,  a  novel  approach  has  been 
adopted  to  isolate  and  to  study  the  variation  of  Reynolds 
stress  components  along  the  flow,  free  from  the 
implications  of  other  modelling  assumptions. 


It  may  be  shown  empirically  that  the  transport 
equations  for  turbulent  kinetic  energy  and  those  for  the 
second  and  third  invariants  of  the  anisotropy  tensor 
contain  information  on  the  distribution  of  energy  among  the 
Reynolds  stress  components.  To  avoid  implications  arising 
from  the  approximate  modelling  of  the  dissipation  equation, 
a  non-d imens i o i na 1  time  r  related  to  the  natural  decay  time 
of  turbulence  energy  is  introduced. 

€  dt  =  dr  ,  with  r  =  — lnt  k/k0)  an<j  dt  =  dx 
k  U 

where  k0  is  the  initial  value  of  the  turbulent  kinetic 
energy.  Then  the  equations  for  the  turbulent  kinetic  energy 


and  the 
expressed 

invariants  of 
as  foil ows : 

the  anisotropy  tensor 

-r 

can  be 

II 

*sj  *- 

k  =>  k  = 

k0  e 

(2  ) 

ai  i 
dr 

=  -2 [ ( c ! -1  )  I  I  - 

X  III) 

(3) 

a  1 1 1 

dr 

=  -3[ (ci-1 ) I I I 

-  X  I  12  ) 

6 

(4) 

where  III  =  b  i  k  b  k  j  b  j  i  ,  the  third  invariant  of  the 
anisotropy  tensor.  By  taking  the  normalised  time  r  from 
experiments,  the  problem  is  transformed  into  the  solution 
of  three  scalar  equations  with  the  normalised  time  as  the 
independent  variable.  Thus  by  the  solution  of  Eqs.  (2)  to 
(4)  the  implications  of  the  modelling  assumptions  of  the 
dissipation  equation  have  been  eliminated  and  the  effects 
of  p r e s s u r e- s t r a i n  term  isolated. 

The  second  invariant,  II,  can  be  interpretated  as  a 
measure  of  the  state  of  anisotropy  of  the  flow.  The  role 
and  physical  meaning  of  the  third  invariant  in  the 
redistribution  of  energy  among  the  components  of  the 
Reynolds  stress  tensor  can  be  obtained  from  the  solution  of 
the  cubic  equation  for  the  principal  Reynolds  stress 
component  s ,  viz. 

b3  -  1/2  I  I  b  -  1/31  I  I  =  0  (5) 

with  b  representing  the  individual  component  of  the 
anisotropy  tensor.  The  solutions  show  that  the  sign  of  the 
third  invariant  remains  constant  during  the  process  of 
return  to  isotropy.  When  the  third  invariant  is  negative, 
the  two  larger  components  are  giving  energy  to  the  smaller, 
consequently  the  smallest  component  tend  to  increase  its 
energy,  while  when  the  third  invariant  is  positive,  the  two 
smaller  components  are  gaining  energy  from  the  largest  and 
the  return  to  isotroy  is  slower  than  the  former  case. 

Analytical  solutions  of  Eqs.  (2)  to  (4)  as  a  function 
of  the  normalised  time  were  compared  with  experimental 
data.  By  neglecting  non-linear  effects  (X=o)  and  making  cj 
a  constant,  Rotta’s  (10)  proposal  is  recovered.  The 
solutions  is  shown  to  be  true  only  for  homogeneous 
axisymmetric  turbulent  flows  or  for  flows  with  small 
departure  from  axisyraraetry.  Similar  solution  also  is 
obtained  when  c^  is  not  a  constant  (X=o).  Therefore  any 
improvement  in  the  modelling  of  the  return  to  isotropy  part 
of  the  pressure-strain  terra,  in  expressing  cj  as  a  function 
of  some  relevant  parameters  (11),  is  still  unsatisfactory 
to  model  non- ax i s ymme t r i c  flows.  Optimisation  of  the  value 
of  the  constant  c i  against  the  data  of  Uberoi  (3) 
(anisotropic  decaying  turbulence  submitted  to  axisymmetric 
plan  strain  with  different  contraction  ratios)  showed  that 
the  value  varies  with  contraction  ratios  as  well  as 
development  of  the  flow.  However  Fig.  1  shows  that  a 
constant  value  with  cj  =  2.0  gives  satisfactory  agreement 
with  the  data  of  Uberoi  (3). 

When  non-linear  terra  in  the  p r e s s u r e- s t r a i n  tensor  is 
considered,  the  axisymmetric  flow  experiment  of  Uberoi  (3) 
can  be  used  to  estimate  the  value  of  the  second  constant  X 
as  a  function  of  that  of  cj .  Figure  2  shows  the  straight 
line  relationship  between  quantity  |X/II0|  and  ( c  j  - 1  )  . 
Employing  different  values  for  c^  and  X  from  Fig.  2,  the 


non-linear  model  predictions  of  Uberoi  (3)  experiment  are 
shown  in  Fig.  1  where  an  improvement  with  respect  to  the 
linear  model  can  be  noted. 

For  non- ax i syrame t r i c  flows,  an  approximate  analytical 
solution  was  obtained  using  perturbation  methods.  The 
solutions  indicate  that  X  and  c)-l  are  important  parameters 
which  influence  the  speed  at  which  the  return  to  isotropy 
is  taking  place.  A  new  parameter  which  measures  the 
initial  degree  of  departure  from  axisymmetry  also  can  be 
defined  as: 

Q  =  1  -  6  I  I I02  I  I0-~  ( 6  ) 

so  that  Q  =  0  for  axisymmetry  and  0-»l  as  the  departure  from 
axisymmetry  increases.  On  an  empirical  basis,  the  speed  of 
the  process  of  equilibration  of  energy  can  be  expressed  by 
the  relationship  between  cj-1  and  the  third  invariant  as 

-  fO  I  I  I  o 

ci  -  1  =  e  (?) 

with  f  being  a  constant  with  a  value  of  10. 

Three  return  to  isotropy  experiments  were  used  to 
compare  with  the  perturbation  and  the  linear  model 
solutions. 

Figure  3  shows  the  comparison  for  the  Tucker  and 
Reynolds  (4)  experiment  where  homogeneous  grid  generated 
turbulence  was  subjected  to  a  uniform  plane  strain.  Good 
agreement  is  obtained  with  the  perturbation  solution  in 
particular  for  the  smallest  component  which  increases  its 
energy  along  the  parallel  duct,  demonstrating  the 
superiority  of  the  non-linear  model. 

Gence  and  Mathieu  (5)  considered  the  return  to 
isotropy  phenomenon  after  the  homogeneous  turbulent  flow 
has  been  subjected  to  two  successive  plane  strains  whose 
principal  axes  are  shifted  by  an  angle  a.  Figure  4  shows 
that  when  ot  =  o  and  a=rt/4,  excellent  agreement  is  obtained 
with  perturbation  solution.  When  oc  =  3tt/4,  the  perturbation 
solution  becomes  less  satisfactory  but  this  may  be 
attributed  to  the  uncertainty  in  applying  the  values  of  cj 
and  X.  which  have  been  found  from  initial  conditions,  to 
the  overall  state  of  the  flow. 

Le  Penven  et  al  (6)  observed  the  evolution  of 
turbulent  kinetic  energy,  second  invariant  and  the  Reynolds 
stress  components,  under  the  influence  of  the  third 
invariant  with  different  signs.  Figure  5  shows  that  when 
the  third  invariant  is  negative,  better  agreement  is 
obtained  with  perturbation  solution  then  with  the  linear 
model,  especially  for  the  smallest  component.  When  the 
third  invariant  is  positive  both  solutions  seem  to  be 
unsatisfactory,  but  again  there  appears  to  be  some 


uncertainty  in 


the  initial  time  of  decay. 

in  conclusion,  a  breakthrough  has  been  achieved  in  the 
understanding  of  the  dynamics  of  energy  transfer  and  the 
non-linear  effects  in  the  return  to  isotropy [  »poce8*  of  Js 
nressure-str ain  term.  The  use  of  non-linear  term  is 
»p'ctVd  to  eliminate  tl>.  nonsttaint  of  axlsymmetry  imposed 
by  linear  modelling  of  Rotta  (10). 

The  work  is  clearly  unfinished  due  to  administrative 
difficulty.  It  was  proposed  to  model  homogeneous 
turbulence  with  mean  strain,  especially  that  of  a  rotating 
turbulence  in  order  to  study  the  functional  dependence  of 
the  constants  and  also  the  non-iinear  effects  in  the  r aP 
term'  or  i  s o t r o p i z a t i o n  of  production  process  in  the 
Pressure  strain  tensor.  Then  the  Reynolds  stress  model 
would  be  employed  to  model  compressible  flows  and  swirling 
flows  of  engineering  relevance. 
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g.  NEW  DISCOVERIES 

The  work  carried  out  can  be  regarded  as  an  important 
breakthrough  in  Reynolds  stress  modelling.  We  are  now  able 
to  model  highly  anisotropic  turbulent  flows  with  a  clearer 
understanding  of  the  non-linear  effects  and  functional 
dependence  of  the  values  of  the  constants. 


h.  ANY  OTHER  STATEMENTS 


See  section  g. 
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